Multiple electron-hole scattering effect on quasiparticle properties 
in a homogeneous electron gas 
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We present a detailed study of a contribution of the T matrix accounting for multiple scattering 
between an electron and a hole to the quasiparticle self-energy. This contribution is considered 
as an additional term to the GW self-energy. The study is based on a variational solution of the 
T-matrix integral equation within a local approximation. A key quantity of such a solution, the 
local electron-hole interaction, is obtained at the small four-momentum transfer limit. Performed by 
making use of this limit form, extensive calculations of quasiparticle properties in the homogeneous 
electron gas over a broad range of electron densities are reported. We carry out an analysis of 
how the T-matrix contribution affects the quasiparticle damping rate, the quasiparticle energy, the 
renormalization constant, and the effective mass enhancement. We find that in comparison with 
the GW approximation the inclusion of the T matrix leads to an essential increase of the damping 
rate, a slight reduction of the GW band narrowing, a decrease of the renormalization constant at 
the Fermi wave vector, and some "weighting" of quasiparticles at the Fermi surface. 

PACS numbers: 71.10.-w 
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I. INTRODUCTION 

It is well known that the extensively used GW approx- 
imation (GWA) to the quasiparticle self-energy is one of 
the most successful methods for describing the quasipar- 
ticle properties in a broad spectrum of materialsiiSi^ 
This approximation employs only the first term in the 
Hedin self-energy diagrammatic expansion^ in the dy- 
namically screened Coulomb interaction (or the so-called 
test charge — test charge screened interaction). However, 
in a number of cases, for example, ferromagnetic tran- 
sition metals," the GWA is less satisfactory. It requires 
the cases to go beyond the GWA (see, e.g., Refs.QHIa)- 

One of the approximations that allows one to relatively 
simply go beyond the GWA is the so-called GWK ap- 
proximation (A stands for the three-point vertex func- 
tion). The latter includes vertex corrections to the GW 
self-energy in the same way as it can be done for the ir- 
reducible polarizability by means of the spin-symmetric 
local-field factor. Thus the GWK approximation, as well 
as the GWA, accounts for charge-density fluctuations 
only. Moreover, this approximation depends weakly on 
the local-field factor and gives results close to those in 
the GWA (see Refs. [l3 and[ll|). Note that it can be 
considered as a GW^-like approximation, where the GW 
formula for the quasiparticle self-energy is used but with 
the electron — test charge screened interaction instead of 
the test charge — test charge one. 

The GW formula for the self-energy can also be de- 
rived in approximations using the exact Ward identity as 
a starting point (see, e.g., Ref . IT^ . In this case, the role 
of the screened interaction W is played by the effective 
electron — electron screened interaction that includes the 
local-field effects by means of the spin-symmetric and 



spin-asymmetric local-field factors and, thereby, takes 
into account the contributi ons of both charge and spin 
fluctuations (see also Refs. Il3ll4ll5il^ and references 
therein) . 

The specific feature of the approaches mentioned above 
is the use of the local-field factors which are defined out- 
side the scope of these approaches and, as a rule, are tab- 
ulated and parametrized by using quantum Monte Carlo 
(QMC) calculations for the homogeneous electron gas. It 
causes certain difficulties in the description of quasiparti- 
cle properties, since these factors do not include real sys- 
tem band structure effects. Thus, it is important to find 
a feasible and all-sufficient scheme for approximate calcu- 
lations of the self-energy which would allow us to preserve 
the advantages of the GWA and at that to give a pos- 
sibility of an inclusion of both charge- and spin-density 
fluctuations, that is crucialiSiiLi^ to obtain agreement 
with experimental data. 

In the preceding paper^i^ hereafter referred to as I, 
we have examined a possibility to go beyond the GWA 
by summing an infinite number of ladder diagrams of 
the Hedin self-energy expansion. To this end, we have 
found a variational solution of the Bethe-Salpeter equa- 
tion (more precisely the ladder approximation to this 
equationSS,) determining the four-point T matrix that de- 
scribes multiple scattering of propagating particles. The 
solution has been obtained within a local approximation 
in the spirit of Refs. and 0. Making use of this so- 
lution, we have proposed a form for the quasiparticle 
self-energy which allows one to take into account charge- 
and spin-density fluctuations without double counting. 
A key quantity of such an approach is the local inter- 
action W . In fact, the sum of the ladder diagrams is 
reduced to this interaction just as vertex corrections to 
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the random phase approximation (RPA) polarizabihty 
can be reduced to the spin-sym metric local-field factor 
(see, 



e.g., Refs. I21l22l23l24l25j) . However, the point is 
that the interaction W is defined by an eightfold integral 
and as such is problematical to be realized in ab initio 
calculations for real systems. 

Thus, the goal of this paper is twofold. Having re- 
stricted ourselves to the multiple electron-hole scattering 
case, first, we find an expression for W suitable for ab 
initio calculations. As we are mainly interested in small 
quasiparticle excitation energy, we consider the small 
four-momentum transfer limit for this interaction^ Sec- 
ond, making use of the obtained limit form of W, we 
analyze what effect the T-matrix contribution being an 
additional term to the GW self energy has on quasipar- 
ticle properties determined by both the self-energy and 
its derivatives with respect to momentum and frequency. 
In this paper, we examine the T-matrix contribution as 
applied to the homogeneous electron gas in the param- 
agnetic state. 

The paper is organized as follows. In Sec. ^ we de- 
rive explicit expressions for the spin-diagonal and spin- 
non-diagonal parts of the local electron-hole interaction 
at the small four- momentum transfer limit. Within the 
model of the homogeneous electron gas in the param- 
agnetic state, by making use of the connection between 
the local interaction and the exchange part of the local- 
field factor, we examine these expressions by comparing 
with the results known from the literature. In Sec. Illll 
we present our main results of extensive calculations car- 
ried out for quasiparticle properties over a broad range of 
electron densities (for values ranging from 2 to 56). On 
the base of these results, we analyze how the T-matrix 
inclusion with the obtained limit form for the local inter- 
action modifies quasiparticle properties evaluated from 
the GW calculations. Finally, the conclusions are given 
in Sec. llVI Unless stated otherwise, atomic units are used 
throughout, i.e., = h = m = 1. 



II. LADDER APPROXIMATION 

In this section, we derive basic expressions defining 
the local electron-hole interaction at the small four- 
momentum transfer limit. Using the relation between 
this interaction and the exchange part of the local-field 
factor, we compare our results for the homogenous elec- 
tron gas (HEG) with those existing in the literature. 



A. r matrix 

The T matrix shown in Fig. ^a) allows one to treat 
the problem of summation of infinite classes of ladder 
diagrams in the diagrammatic expansion of both the ir- 
reducible polarizabihty [Fig. ^b)] and the quasiparticle 
self-energy [Fig. ^c)]. In momentum space, within the 
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FIG. 1: A diagrammatic representation of the T matrix (a), 
accounting for multiple electron-hole scattering, the consid- 
ered ladder diagrams for the irreducible polarizabihty P (b), 
and the T-matrix contribution S"^ to the self-energy (c). (d) 
Feynman diagrams for Kc„i[Q) and M„„i{Q). The wiggly 
lines signify the dynamically screened Coulomb interaction 
W . The solid lines with arrows represent the Green function 
G. The T matrix is shown by the shaded square. 



local approximation the T matrix accounting for multiple 
electron-hole (e — h) scattering has the formiS 



'{Q) 



l-Waa,{Q)K,^,{Q) 



(1) 



Here and in the following we use the four-momentum 
variables Q, k, or p as a shorthand for (Q,r2), (k, cj), 
or (p,S7), respectively, tr labels the spin, and a — ± 
corresponds to spin-up, and spin-down, respectively. The 
e — h propagator K^^i is given by^^ 



-GAQ+p)Ga'{p), 



(2) 
(3) 



where Go- is the Green function. The local electron-hole 
interaction Waa' is related with the dynamically screened 
Coulomb interaction W by the equation [see Fig. ^d)] 

W^a'iQ) = [Kaa'{Q)VM,,,{Q)[K„„,{Q)]-\ (4) 
where 

M„„,{Q) = J dkdp 

{k)W{k-p)K^^,M- (5) 

Note that the T matrix depends only on the four- 
momentum transfer along the electron-hole channel. 

In spite of this simple form, a calculation of the T ma- 
trix for the real system is still difficult due to the eight- 
fold integration in Eq. defining the local interaction. 
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In this paper, in order to obtain a suitable expression Ko- -a,Qik) we can use the ratio 



for Waa' , we consider the small four-momentum transfer 
limit (i.e., we set f2 = and then take the limit Q ^ 0). 
It is expected to be a reasonable approximation because 
we are mainly interested in excitations in the vicinity of 
the Fermi energy where quasiparticles are well defined. 



B. Local interaction 

First of all we examine the spin-diagonal part of the lo- 
cal interaction. Omitting argumentations which one can 
find in Ref . '27", in the calculation of poles contribution to 
the integral lO we replace Ko-o-,g(fc) by CaaS{uj)S[ec{'k)]. 
Unless stated otherwise, the quasiparticle energy (k) is 
measured from the Fermi energy ep- The coefficient Caa 
is determined by the relation |(2Jl. Thus, Ha<T,Qik) can be 
written as 



..Q(fc) 



1 



:K^^{Q)5{uj)S[e^{k)]+^{k), (6) 



where (p{k) is the regular part of Kc,cr,Q{k) and Np is 
the density of states of spin a per unit volume at the 
Fermi surfacei^ Substituting Eq. JHl into Eq. © and 
neglecting2S terms with ip, we obtain that at the consid- 
ered limit the local interaction has the form 



or 



f dkdp6[e4k)]W{k - p, 0)S[e^{p)] 



-^W^(k-p,0)^ 



[Sn^N^VJs. K(k)|- - "'"'WApW 



(7) 



(8) 



where Sp is the Fermi surface for spin cr, and Vo-(k) = 
Vkeo-(k) is the quasiparticle velocity at this surface. Now 
we have the fourfold integration instead of the eightfold 
one in Eq. (0). 

For the HEG in the spin-polarized state, Eq. ijSJl can 
be rewritten as 
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(47r) 



df7kdf2pTV(fc|;|k-p|,0), 



(9) 



where fJk and Qp are spatial angles, k and p are unit 
vectors. The Fermi wave vector kp for spin a is related to 
that in the paramagnetic state kp by hp = kpi^ — crQ^^^, 
where the relative spin polarization ^ = \n-^ — n-\/n (jia 
is the spin a electron density, ri = ri+ + n_ being the 
HEG electron density) defines the exchange splitting of 

the band as 2 A = -^[(1 + C)^^^ - (1 - 0^^^]-^ Note 
that considering the HEG, here and in the following, in 
all expressions we use the noninteracting energy (k) — 
k^/2 — kp"^ /2 in place of ecr(k), but for the quasiparticle 
energy in this case the notation Ek is used (see Sec lIIIBl . 

As to the spin-non-diagonal part of the local inter- 
action, instead of (5[eo-(k)] in the expression replacing 



rtf [e^(k)] -nf [e-^(k)] 
^-^^''^^ .-.(k)-e.(k) ' (^0) 

where np is the Fermi distribution function. In this case, 
the coefficient Ca -a is equal to Ka -a{Q)l ^■n'^Ka -o- (0)] . 
As a result, the spin- non-diagonal local interaction is 
given by the form 



J dkdp7^_^(k)iy(k-p,0)7^_^(p) 
[87r3if,_,(0)]2 



(11) 



For the HEG in the spin-polarized state, we can rewrite 
Eq. Ijllfl in the following way 



1 



(47r)2 

k J-. 



1 



[2C4/3]^ 



(12) 



Ikrdlkl 



k1 



|ppd|p|Ty(k-p,0). 



Note that the ratio ja-a tends to the (5-function at 
the C — > limit (the paramagnetic state), and the 
spin-diagonal and spin-non-diagonal parts become equal. 
Since in this paper we are interested in paramagnetic 
systems, the local interaction with C remains to be 
examined elsewhere. 



C. Local-field factor 

In paper I, by examining the irreducible polarizability 
exchange diagrams [Fig. H^b)], we have shown that for 
paramagnetic systems the local electron-hole interaction 
W can be identified with the local-field factor G{Q) (or 
more precisely with its exchange part) 



GiQ) = W{Q)/2v,{Q), 



(13) 



where W{Q) = ^J2a^("^iQ)^ ^^'^ ^ciQ) is the bare 
Coulomb interaction. This exchange part is related to 
the first-order exchange diagram in the irreducible polar- 
izability diagrammatic expansion [see Eqs. Q), (0) and 
Figs.n^b) andn(d)]. Such a relation was derived and ex- 
amined by m any authors wit hin var ious approaches (see, 
e.g., Refs. l2lll24i25l and'3l"3 ll33l34|) . Thus we can verify 
the approximation done for W in the previous subsection. 

Actually, at the considered limit {Q 0), the local- 
field factor of the HEG has the form 



Q 



(14) 



where A — ^^-^ in our case. This relation allows us to 
carry out, by the example of the HEG in the paramag- 
netic state, a comparison of the coefficient A obtained 
with the help of Eq. (O with that known from the liter- 
ature. 
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FIG. 2: The coefficient A versus the electron density param- 
eter ra. Filled squares and circles represent the results of 
Refs. and respectively. Stars with the error bar de- 
pict A obtained from calculations of Qa for rs = 2 shown in 
Fig. 3 of Ref.m The notations H, CA/PZ, TF, L, and W 
signify, respectively, the Hubbard approximation, the results 
of Monte Carlo calculations of Ref. (ra > 1), the Thomas- 
Fermi model [Eq. 1151 1. calculations by means of Eq. ® with 
the use of the Lindhard function, and, finally, A obtained 
from Eq. 1161 . Open triangles and diamonds depict the "self- 
consistent" calculations (see the text) for the TF and W cases, 
respectively. 

The simplest way to analytically perform the integra- 
tion in Eq. is to use the screened Coulomb interac- 
tion VF(k,0) = 477 [|k|2 of the Thomas-Fermi 

model (see, e.g., Ref. [s^. Here qtf = {^^^Y^^kp is 
the Thomas-Fermi wave vector, the electron density pa- 
rameter rs is given by ^{aorg)^ — l/n, qq being the 

1 /3 

Bohr radius, ce — {■^) , a-nd argkp — 1. As a result, 
for A we obtain the following dependence on the electron 
density parameter: 

ATF{rs)^lln(l + —) . (15) 

In Fig. 121 we show Atf as a function of rs as well as 
the coefficient Al (the "L" curve) calculated by means 
of Eqs. ijni, CSli and ifTH) with the use of the Lindhard 
dielectric function. It follows from the figure that Atf 
and Al are very close to each other, especially at 0. 
Moreover, both these results are in good agreement with 
the coefficient A obtained in Ref. |2J (in the figure it is 
shown for Al, Li, Na, and K), where the same class of 
ladder diagrams in the irreducible polarizability diagram- 
matic expansion was considered. 

Also, as shown in Fig. [2 our results are in accord with 
calculations performed in Ref. I21I for the local-field factor 
Qa [sec Eq. (36) in the cited paper] which is formally the 
same as that of Eq. H13|) , except for including self-energy 
effects and the fact that Qa is given on the imaginary- 
frequency axis. The shown error bar is originated from 
the proc edure of evaluating of A: inspecting Fig. 3 of 
Ref. I21I we have found A from a set of values of Qa 



evaluated at various momenta {OAkp, O.Gkp, O.Skp, and 
IMf). 

Due to the integration over the Fermi surface in 
Eq. lO, at the low-density limit (ts — > 00), when the 
ratio qxF/kp is large, one can neglect momentum depen- 
dence of the screened Coulomb interaction. This leads to 
W = W{Q, 0) and, as a consequence, 

Aw{rs) = (16) 
aars 

It is appropriate to mention here that the approximation 
made in Ref. UtI corresponds to the local interaction given 
by W{Q) = VF(Q,0). This means that the local-field 
factor of Eq. H13|) can be represented as 

For the HEG, the approximation H17|l gives a Hubbard- 
like form for Q and at the small four-momentum transfer 
limit it leads to A — Aw . Such an approximation can be 
valid at the low-density limit, whereas for Ts < 2.5 it es- 
sentially overestimates the exchange diagrams contribu- 
tion, exceeding even the accurate results of Monte Carlo 
calculations^ which include both exchange and correla- 
tion effects (see Fig. 

It is worth pointing out that the common property of 
Atf, Al, and Aw is the divergence at Vs — > 0."^^ In this 
connection, it makes sense to compare our results with 
the exchange local-field factor of Ref. 22 (filled squares 
in Fig. EJ. A class of diagrams considered there com- 
prises ladder diagrams based on more complex electron- 
hole interaction containing electron-electron and hole- 
hole multiple-scattering events. As a result, the coeffi- 
cient A of Ref. ^ (evaluated from the local-field factor at 
kp shown in Fig. 7 of the quoted paper) has no divergence 
in its dependence on rs , tending to the Hartree-Fock pre- 
diction 1/4 at rs 0. Nevertheless, by inspecting Fig.|21 
we can infer that for the entire metallic density range (r^ 
from 2 to 6) even the simplest approximation At f yields 
results consistent with those known from the literature. 
This fact allows us to expect that the local electron-hole 
interaction {T))--® will give reasonable results on the self- 
energy ladder diagrams treatment. 

At last, we would like to note that in order to more pre- 
cisely evaluate the exchange part of Q from Eqs. © and 
(|13|l , one can include the corresponding local-field correc- 
tions into W and perform a kind of "self-consistent" (sc) 
procedure. '^^ As shown in Fig. |21 (open triangles and dia- 
monds), such a procedure modifies Atf slightly, whereas 
in the case of Aw the changes are more significant and 
the resultant Aw-.sc = ^Aw allows one to reasonably es- 
timate the considered diagrams contribution in the entire 
metallic density range (^Moreover, the use of this prefac- 
tor 2/3 together with the approximation H17|l in ab initio 
calculations^" brings theory and experiment to a better 
agreement. 
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III. QUASIPARTICLE PROPERTIES 

In this section, we address the question of how 
the electron-hole multiple scattering affects quasiparticle 
properties which are determined by both the self-energy 
and its derivatives with respect to momentum and fre- 
quency. We study the HEG in the paramagnetic state at 
different values of ranging from 2 to 56. This study 
includes the metallic density range (2 < Ts < 6) as well 
as the range of large Vs ^ 48 for which the effective mass 
"divergence" was predicted in Ref. li^ . 

A. Quasiparticle self-energy 

In paper I, we have shown that the T matrix al- 
lows one to go beyond the GWA by summing an infinite 
number of the electron-hole ladder diagrams shown in 
Fig. ^c). In this case, the quasiparticle self-energy can 
be expressed as T,^ = -|- Ej, where the GW term is 
well known to be 

S^^b) = / dkGAk)W{p - k), (18) 

and the T-matrix contribution is given by 

^lip) = -(^E / dkG,,{k)%,,{p^k) (19) 

with %a'{k) = taa'{k)[K^^,{k)Waa'{k)f- It is impor- 
tant that the T-matrix contribution (|19|) has a GM^-like 
form that simplifies calculations of Ej. 

Due to the correspondence between multiple-scattering 
events and a spin fluctuation (see, e.g., Refs. and 
[T7..43..45-46..47.48) . this approach to the self-energy in- 
cludes the contributions of both charge and spin fluctua- 
tions. As a result, in such an approach, having retained 
all the advantages of the GW approximation, we have 
a possibility of describing quasiparticle properties more 
comprehensively than it can be done in the GWA. Note 
that as in the case of the irreducible polarizability lad- 
der diagrams, the local interaction W^a' is an object of 
principal concern here because T^o-' depends significantly 
on the form for this interaction and, consequently, care 
must be taken by choosing an approximation to it. 

In this work, we mainly focus our attention on the ef- 
fect of the inclusion of the self-energy ladder diagrams 
on the quasiparticle properties in the case of the HEG in 
the paramagnetic state. For simplicity, we evaluate the 
GW term with the noninteracting Green function and the 
RPA screened interaction (the so-called G^Wq approxi- 
mation which gives a better description of the quasiparti- 
cle properties than the fully self-consistent GW approx- 
imation itself, except for the total energy^ii^). Below, 
unless stated otherwise, this GqW^ approximation is re- 
ferred to as the GWA. To evaluate the T-matrix contri- 
bution, we also use the noninteracting Green function. 



the RPA electron-hole propagator, and the local interac- 
tion elaborated in the previous section (corresponding to 
Al). 

For the computational purposes, it is convenient to 
split up the GW term into the energy-independent 
Hartree-Fock (exchange) part E^-'^(k) and the corre- 
lation part S^^(k,u;) defined through Eq. (HHJ) by 
the induced potential Wi = W — Vc instead of the 
screened interaction W. Within the spectral function 
representation,^ having found the imaginary part of E^*^ 
we obtain its real part, ReSjf^(k, w), from the Hilbert 
transform by using the principal value integration. As a 
result, the real part of the self energy in the GWA is equal 
to ReE^^(k,tj) = S^^(k) + ReJ:^^ {k, lu) . Regard- 
ing the T-matrix contribution, we first evaluate ImS^ of 
Eq. H19() and then perform the Hilbert transform. 

In order to get some idea of the quantity and behav- 
ior of the T-matrix contribution in comparison with the 
GW term, in Fig. O we plot the correlation part of the 
self-energy and the T-matrix contribution E"^ as a 

function of the momentum and energy for = 4. As fol- 
lows from the figure, due to plasmon singularities of the 
inverse dielectric function, the GW term seems substan- 
tially larger than the T-matrix contribution.^" However, 
inspecting Figs.|3fc) and|2Id), one can see that the 
surface becomes notably corrugated in the vicinity of the 
uj = energy where T,^^ behaves smoothly and contains 
only the electron-hole contribution (the plasmon decay 
channel is not open yet in this energy range) which is 
comparatively not large. As a consequence, the T-matrix 
inclusion is expected to appreciably influence the quasi- 
particle properties (especially on those determined by the 
self-energy derivatives) in this energy region. 

As is evident from Fig.|2Ic), at some energies co < the 
spectral function of the T-matrix contribution becomes 
negative. However, it does not lead to wrong analyti- 
cal properties of the self-energy, because the contribution 
(|19|l has such a form as an additional term beyond the 
GW one and, consequently, should be considered in sum 
with the latter. This sum is always non-negative. 

Now we consider the on-shell imaginary part of the 
self-energy which gives the quasiparticle damping rate. 
In Fig. ^ we show this quantity as a function of mo- 
mentum jkl/Zci? for — 2 and 4. First of all, we would 
like to note that despite the fact that the spectral func- 
tion of S^(k,w) is negative within some k — w domain 
the on-shell imaginary part of the T-matrix contribu- 
tion has the same sign as the GW term at any value 
of momentum. This means that the T-matrix inclu- 
sion leads to an increase of the quasiparticle damping 
rate with respect to that in the GWA. In order to as- 
sess the quantity of this increase, we examine the ratio 
ImE/ImS'^^ = 1 + ImS'^/ImE^^ (see inset of Fig-EJ. 
According to the calculations performed, this ratio as a 
function of momentum has the largest value at |k| kp 
and falls down when |k| moves away from kp. Fur- 
thermore, the lower the electron density that the HEG 
possesses, the larger the ratio we have. For example. 
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FIG. 3: The spectral function — sgn(aj)ImE;, (a) and the real part ReE;, (b) of the correlation part of the self-energy 
calculated within the GWA. Also the spectral function — sgn(aj)ImE"^ (c) and the real part ReS"^ (d) of the T-matrix contribution 
to the self-energy. The GWA self-energy and the T-matrix contribution are plotted as functions of energy Loj^F and momentum 
|k|/fcF for ~ 4. At the top of each of these figures, the contour plot (50 levels) is shown. 



ImE/ImE'^^ evaluated at |k| = 1.05kp is equal to 1.33 
for Ts — 2 and 1.72 for = 4. 

As compared with the generalized GW self-energy 
evaluated in Ref. l4ll for |k| < kp we have similar results 
which are very close quantitatively. However, it is not so 
for |k| > hp especially in the region where a quasiparticlc 
can decay into plasmons. Under this region the on-shell 
imaginary part of the generalized GW self-energy very 
quickly becomes smaller than that in the GWA and then 
demonstrates a partial "suppression" of the plasmon de- 
cay channel in comparison with the GW approximation 
(see also Ref. lsij) . 

Next, to answer the question of how the value of the 
ratio ImE/ImE*^^ depends on the chosen form of the 
approximation for the local electron-hole interaction, we 
have calculated ImE'^ with W defined by Eq. H13|) with 
the local- field factor (|17|l for rg — 2. As follows from the 



inset of Fig. 01 the ratio becomes essentially larger than 
that presented by dashed line and runs up to 2.53 at |k| = 
l.Obkp- This means that the use of such an interaction 
W instead of Eq. ^ leads to an increase of the ratio 
by a factor of 1.90 at this momentum. Thus, keeping 
in mind that at the small four-momentum transfer limit 
this local interaction corresponds to the coefficient A\y 
of Eq. (|16(l . we suppose that at least for aluminum the 
use of this approximation to W results in a large value of 
the T matrix that in turn rises the quasiparticlc damping 
rate too high. In connection with the coefficients shown 
in Fig. 121 we would hke to note that, e.g., for = 2 
the local-field factor of Ref. |2^ used in paper I gives the 
increase of the ratio by a factor of 1.24. 

Regarding the on-shell real part of the self-energy (see 
Fig. [SJ, we would like to emphasize that again the T- 
matrix contribution has the same sign as the GW term. 
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FIG. 4: The absolute value of the imaginary part of the 
self-energy evaluated at lj = e(k) (the GW term and the 
T-matrix contribution) as a function of momentum |k|/fcF 
for = 2 (dashed line) and 4 (solid line). Inset: the ratio 
ImE/ImE'^^ = 1 + ImE^/ImE"^^ as a function of momen- 
tum |k|/fcF for rs = 2 and 4. The dash-dotted line represents 
this ratio evaluated with the local-field factor of Eq. I|17|l for 
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FIG. 5: The real part of the self-energy ReE[|k|, e(k)] (the 
GW term and the T-matrix contribution) as a function of 
momentum |k|/fcF for rs = 2 (dashed line) and rs = 4 (solid 
line) . 



Owing to the "corrugations" discussed in connection with 
Fig. 121 in the vicinity oi kp ReS"'" shows fast variation 
which for 0.7 kp < |k| < l.Skp can be fitted rather well 
by a sinehke function. As foUows from Fig. |S1 the T- 
matrix contribution represents itself as a very small quan- 
tity within all the considered k-domain and therefore can 
alter the on-shell real part of the self-energy evaluated in 
the GWA to only a small extent. 

In Figs. El and [7| we show how the GW term and the 
T-matrix contribution depend on beyond the metallic 
density range. One can see that the on-shell ImE'-^^ 
as a function of shows relatively small changes for 
O.Qkp < |k| < 1.5kp and already at Vg ^ 24 gets 
some "saturation" , whereupon properties of the e — h 
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FIG. 6: The on-shell real and imaginary parts of the self- 
energy calculated in the GWA as functions of momentum 
jk|/fcF for Vs values ranging from 4 to 56. 
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FIG. 7: The on-shell real and imaginary parts of the T-matrix 
contribution as functions of momentum |k|/fcF for rs values 
ranging from 4 to 56. 



decay channel remain practically unchanged within the 
mentioned k interval. For |k| > l.bkp, especially in 
the plasmon emission region, the quasiparticle damping 
rate demonstrates a continual increase with increasing 
rs- Owing to this and the rg-dependence of E^^, the 
on-shell real part of the GW self-energy monotonically 
decreases as a function of Vg- As to the T-matrix contri- 
bution, after ~ 40 the changes of the imaginary and 
real parts become insignificant for any momentumcS 

In Fig. IHl we plot the ratio ImE/ImE'-^^ as a func- 
tion of Ts evaluated for different momenta. In fact this 
figure shows how the inclusion of the self-energy ladder 
diagrams affects the quasiparticle damping rate. The 
biggest contribution of these diagrams is observed in the 
vicinity of kp and can exceed the GWA prediction, e.g., 
by a factor of 3.6 at Vs = 56. 
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FIG. 8: The ratio ImE/ImE'^^ = 1 + ImE^/ImE°^ as a 
function of ra at |k| = O.bkp (dashed Hne), l.Obkp (soUd 
line), and l.bkp (dotted line). 



B. Quasipartical energy, renormalization constant, 
and effective mass enhancement 



The quasiparticle excitation energy within the first- 
order perturbation theory is given by 

Ek ^ e(k) +ReI][k,e(k)] -ReE[A:i.,0], (20) 

where the real part of the self-energy in the on-shell ap- 
proximation is added to the single-particle energy. In 
Fig. O we show the calculated Ek for Vg = 2 and 4. 
Due to the small value of the on-shell real part of the T- 
matrix contribution, the quasiparticle dispersion is only 
slightly changed by the inclusion of the self-energy ladder 
diagrams. Inspecting the figure, one can see that this in- 
clusion results in a light increase of the bandwidth with 
respect to the GW one. This increase is similar to that 
evaluated in Ref.T^ within the GWA approximation with 
the inclusion of the same vertex function in the screened 
interaction and the numerator of the self-energy. How- 
ever this similarity is observed only under Ts ~ 4. For 
''s si 4, the GWA approximation yields band narrowing 
greater than the GWA does. 

As was anticipated, we have a more profound effect 
of the T-matrix inclusion on the quasiparticle properties 
determined by the derivatives of the self-energy. One of 
such properties is the effective mass enhancement which 
in the on-shell approximation is known to be given by 



m*(k) 



m dEk 



k dk 



(21) 



The inset of Fig.|51represents the inverse value of this dis- 
persing quasiparticle effective mass for Vs ~ 2 and 4. A 
close examination of m/m* as a function of momentum 
provides interesting insights. First, owing to the sinelike 
behavior of the on-shell ReS-^ in the vicinity of fc^ (see 
Fig. ISJ , the T-matrix inclusion gives alternating contri- 
bution to the inverse effective mass. As a result, m/m* 
becomes smaller at |k| ~ kp and greater away from it. 



,^ 0.0 




FIG. 9: The quasiparticle energy Ek, Eq. H2()|l . calculated 
without {GW, solid line) and with (GW&T, dotted line) the 
T-matrix contribution as a function of momentum |k|/fcF for 
Va = 2 and = 4. Inset: the inverse effective mass enhance- 
ment m/m* as a function of |k|/fcF for rs = 2 and 4. 
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FIG. 10: The real and imaginary parts of the renormalization 
constant Zk calculated without (GW) and with (GW&cT) the 
T-matrix contribution as a function of momentum |kj/fcF for 
= 2 and 4. 



At some momenta the T-matrix contribution to the ef- 
fective mass is equal to zero. Second, for Ts = 2 multiple 
electron-hole scattering modifies the effective mass to an 
extent that quasiparticles become "heavier" than in the 
noninteracting system, whereas the GWA predicts the 
reverse. Third, as can be seen from Fig. |S1 for = 4 
the absolute value of the derivative of the on-shell ReS^ 
with respect to momentum should be larger than that for 
= 2. This entails the larger alteration of m/m*. Thus, 
focusing our attention on the effective mass behavior in 
the vicinity of fc^?, we can infer that the multiple electron- 
hole scattering leads to a "weighting" of quasiparticles. 

Note that Eq. (|21|l is a valid approximation to the ef- 
fective mass enhancement at small values of Vg ^ In order 
to estimate the quasiparticle effective mass at kp more 
precisely, especially for large , we use the formally exact 



TABLE I; The effective mass enhancement m* /m and the 
renormahzation constant Zf at the Fermi wave vector kp 
calculated without {GW) and with (GW&T) the T-matrix 
contribution. In each of these cases, two values of the mass 
enhancement are presented: the first one is obtained from the 
exact formula, Eq. 1221 . the second one (in parentheses) is 
calculated within the on-shell approximation, Eq. I21II . 



' I ' I ' I ' I .| 'q I ' I ' I ' I ' I ' I ' I ' I ' I ' 
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Zf 




GW GWkT 


GW 


GWkT 


2 


0.99(0.99) 1.01 (1.01) 


0.77 


0.75 


4 


1.05 (1.07) 1.08 (1.12) 


0.64 


0.62 



equatio] 



,35 



m 

TO 



1 + 



dk 



(22) 



oj— 0,|k|— /cii 



Here Zp is the renormlization constant Zk evaluated at 
the Fermi wave vector. In its turn, the renormalization 
constant Zk that gives the spectral weight of the quasi- 
particle is defined for a four-momentum [k, e(k)] as2^ 



Zk 



as(k,tj) 



duj 



.=e(k) 



(23) 



In Fig. IIUI we plot the renormalization constant as 
a function of momentum for — 2 and 4. It fol- 
lows from the figure that changes of the real part of 
Zk caused by the T-matrix inclusion occur mainly up 
to |k| ~ LS/cF- The imaginary part of Zk accounting for 
multiple electron-hole scattering differs weakly from that 
in the GWA. At |k| = kp, the renormalization constant 
is real. Due to the fact that the frequency derivative 

is negative as well as in the case of the 

GW term (see Fig. EJ, Zp becomes smaller with the T- 
matrix inclusion. Together with the positive momentum 
derivative of ReS"^ it affects the effective mass of Eq. H22|) 
and moves the latter away from its value obtained in the 
GWA. Table m contains our detailed results on the effec- 
tive mass enhancement and the renormalization constant 
for Ts — 2 and 4. 

In Fig. ^2 we show the calculated dependence of Zp 
and m/m* on and compare it with that known from 
the literature. As follows from Fig. Ilir a). the inclu- 
sion of the T matrix reduces the renormalization con- 
stant Zp. Thereby the T-matrix contribution shifts the 
Zp-cuTve towards the so-called GZ (Gori-Giorgi-Ziesche) 
parametrization of Ref. Is^ This parametrization is in 
good agreement with calculations of Ref. |^ performed 
by making use of the cffective-potential-expansion (EPX) 
method which in contrast to the Hedin expansion^ is 
formulated in terms of the static screened interaction 
VF(q, 0). The GZ parametrization is also compatible with 
QMC data found in Ref. |0 for the HEG momentum dis- 
tribution. The fact that the renormalization constant Zp 
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FIG. 11: (a) The renormalization constant Zf at the Fermi 
wave vector obtained without {GW , solid line) and with 
{GWk.T, dashed line) the T-matrrx contribution as a func- 
tion of Ts- The dotted line represents the so-called GZ 
parametrization (Ref. Wk) which is valid in the density range 
Ts < 12. The inset shows the obtained results in compari- 
son with those known from the literature, (b) The inverse 
effective mass enhancement m/m* at the Fermi wave vector 
calculated both in the on-shell approximation, Eq. H21|l . la- 
beled as "on-shell" , and by making use of the formally exact 
Dyson equation 12211 . labeled as "off-shell". The notations 
OB, NI, GV, RS, and HB sign ify the correspondent values 
taken from Refs. Is^. l53l l4lL l56l andl4gl respectively. 



determined by Monte Carlo^ (labeled as OB in the fig- 
ure) is noticeably larger than the parametrization can be 
explained by the difference in procedures of finding the 
momentum distribution discontinuity at kp . In the QMC 
calculations there is a finite distance between momenta 
(including the closest to kp) for which the momentum 
distribution is calculated (see, e.g., Fig. 12 in Ref. 'sF), 
whereas the GZ parametrization allows one to find Zp 
from |k| kp ±Q limits. 

Note that the self-consistent schemes of Refs. 49 (HB), 
tsS (NI), and Is^ (RS), exhibit an increase of the quasi- 
particle renormalization constant at the Fermi surface 
with respect to Zp"^ shown in Fig. Illf a). It can be 
understood within the framework of the detailed analy- 
sis carried out in Ref. where the self-consistency be- 
tween the Green function and the self-energy within the 
GW approximation has been achieved. As was shown^ 
the self-consistency procedure leads to the strongly sup- 
pressed plasmon peaks in the imaginary part of the self- 
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energy. Owing to the Hilbert transform, it entails an 
essential smoothing of the corresponding sharp struc- 
tures of the self-energy real part. As a consequence, 
the frequency derivative of the latter, which is nega- 
tive, becomes smaller by absolute value, that, in its turn, 
through Eq. ^'2'.^ leads to an increase of the renormaliza- 
tion constant. 

As regards the effective mass enhancement, it is evi- 
dent from Fig. lllf b) that the on-shell approximation H21|) 
underestimates m/m* in comparison with that obtained 
from Eq. (|22|l . At large values of rg it comes into particu- 
lar prominence. The on-shell GW effective mass exhibits 
a divergence at rs ~ 50 (in Ref. at rs ~ 48). The 
T-matrix inclusion results in a considerable increase of 
the on-shell effective mass and leads to the divergence 
at Tg ^ 33. The effect of taking into account multiple 
electron-hole scattering on the effective mass calculated 
by making use of the formally exact equation (|22|l is con- 
sistent with the influence of the T-matrix inclusion on 
the renormalization constant. In this case the quasipar- 
ticle mass demonstrates a relatively weak dependence 
without any divergence up to the largest Vs considered in 
the paper. 

Regarding the dependence of the quasiparticle effec- 
tive mass, one can note that the effective mass is a more 
controversial quantity than the renormalization constant. 
Actually, comparing our m/m*(r.^) with that of calcu- 
lations of Refs. ii ("RS") andlHl ("NI"), we find that 
contrary to our results the self-consistent schemes pre- 
dict a monotonic increase of the inverse effective mass 
as a function of Vg. Such an increase is important to 
imply a bandwidth widening at metallic densities that 
disagrees with the experimental findings. Thus together 
with Zp slightly "overestimated" in comparison with the 
GWA [see inset of Fig. IllT a)] these schemes yield the 
effective mass smaller than that in the noninteracting 
system, whereas the GWA gives the reverse. According 
to Eq. H22|l . it means that in the RS and NI cases the 
momentum derivative of the self-energy is larger than 
that evaluated in the GWA. Seemingly, it is caused by a 
large contribution of the Hartree-Fock part which is the 
self-energy for the noninteracting systems providing the 
effective mass equal to zero at the Fermi surface. This 
contribution cannot be canceled by the correlation part 
of the self-energy in full measure, as it occurs in the non- 
self-consistent GWA (see Ref. l49|) . As a result, this leads 
to an increase of the m/m* ratio. To all appearances, 
at least for Ts < 4 we have a similar situation in the 
generalized GW self-energy calculations of Ref. 0. 



IV. CONCLUSIONS 

In conclusion, we have presented a detailed study of the 
effect of multiple electron-hole scattering on quasiparti- 
cle properties determined by both the self-energy and its 
derivatives with respect to momentum and frequency. To 
take into account multiple scattering between an electron 



and a hole in calculations of the self-energy S, a varia- 
tional solution of the corresponding Bethe-Salpeter equa- 
tion obtained in our preceding pape^i^ has been used. 
This solution representing the T matrix within a local ap- 
proximation allows one to sum the self-energy ladder dia- 
grams. To preserve the advantages of the GWA, we have 
considered the sum of these diagrams as an additional 
term (the T-matrix contribution E-^) to the GW self- 
energy E'^'^ . In this approach all weight of the problem 
is transferred^to a form chosen for the local electron-hole 
interaction W that appears in the definition of the T ma- 
trix. By examining the irreducible polarizability ladder 
diagrams, one can identify this interaction with the ex- 
change part of the many-body local-field factor. Consid- 
ering this local interaction at the small four-momentum 
transfer limit, we have arrived at the expression which 
gives the results for the exchange local-field factor in ac- 
cordance with those known from the literature. 

Using the obtained form of W ^ we have carried out 
extensive calculations of both the Yj'-^^ and Yj^ terms 
and such quasiparticle properties as the damping rate, 
the quasiparticle energy, the renormalization constant, 
and the effective mass enhancement over a broad range 
of electron densities in the homogeneous electron gas. 
The calculations have shown that the T-matrix inclusion 
leads to an increase of the quasiparticle damping rate es- 
pecially in the vicinity of kp. This increase depends on 
Ts and can exceed the GW prediction by a factor of 1.8 
for Ts = 12 and 3.6 for rg = 56. Regarding the question 
of how a form chosen for the local interaction affects the 
quasiparticle damping rate, we have found that the latter 
is a very form-sensitive quantity, and consequently it is 
easy to over(under)estimate the ladder diagrams contri- 
bution. We have also revealed that due to small values of 
the on-shell real part of H'^ the T-matrix inclusion modi- 
fies slightly the quasiparticle dispersion reducing the GW 
band narrowing. 

We have found that the T-matrix contribution can no- 
tably affect the renormalization constant and the effec- 
tive mass enhancement. Examining the renormalization 
constant Zp as a function of Vg at the Fermi surface, we 
have ascertained that in comparison with the GW values 
the T-matrix inclusion reduces Zp. As a result, the latter 
becomes closer to the renormalization constant given by 
the GZ parametrization^ compatible with QMC calcu- 
lations of the momentum distribution. A close analysis 
of the quasiparticle effective mass m* /m evaluated both 
in the on-shell approximation and in the formally exact 
Dyson scheme (the off-shell approximation) has shown 
that the on-shell effective mass depends strongly on Vg 
and has a divergence'*^ which is shifted by the T-matrix 
inclusion from Vg ^ 50 (in the GWA) to Vg ^ 33. The off- 
shell m* /m exhibits relatively weak rg dependence and 
does not diverge up to the largest Vg considered. In this 
case, the T-matrix contribution leads only to a slight in- 
crease of "the weight of quasiparticles" . 
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